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I. INTRODUCTION 



The quantum dynamics of systems that are classically chaotic has been a subject of considerable interest for nearly 
three decades. Chaos, usually defined as the exponential sensitivity of phase space trajectories to slight changes 
of the initial conditions, resists direct translation to the natural Hilbert space setting of quantum mechanics since 
there do not exist well-defined trajectories within the conventional interpretation of quantum mechanics Hence, 
many studies on quantum chaos in the literature have been dedicated to the relatively modest study of quantum 
manifestations of classical chaos |2J,|3|. Actual quantum-classical correspondence (QCC) in classical chaotic systems 
is still one of the outstanding issues in quantum theory. 

Considerations of quantum-classical correspondence that rely upon the standard formulations of quantum and 
classical mechanics are at a great disadvantage. Specifically, classical mechanics does not admit wavefunctions and 
quantum mechanics does not admit trajectories. Hence one is faced with trying to connect two theories which do not 
even have the same basic rudimentary elements. 

A far more enlightening approach is to utilize the Hilbert space formulation of both quantum £| and classical [f| 
mechanics where the central element is the density operator in both formulations. By choosing to deal with the phase 
space representation of quantum mechanics one can then compare the classical and quantum dynamics of distributions 
in phase space. Indeed, all the associated machinery of commutation rules, eigenvalues, eigenstates, etc. can be used 
to formally analyze quantum-classical correspondence of both integral and chaotic systems Q . 

Reco gniz ing that phase space distributions allow for a common view of classical and quantum mechanics (see also 
uH&OIUIl) suggests that we consider the nature of chaos as it manifests itself in phase space distributions, as distinct 
from classical phase space trajectories. To this end we ^lj, and others [T^.[T^|. have developed a criteria for chaos in 
terms of the structure of phase space distribution functions. 

In this paper, the distribution function approach to QCC in classically chaotic systems is carefully examined and 
substantially extended. The extension is from one-dimensional chaotic systems in previous work Il2| to two- 
dimensional conservative systems, from the special case of uniform dynamical instabilities to the general case of 
nonuniform stretching and contraction rates fluctuating with time and phase space location. That is, armed with a 
measure of the structure of phase space distribution functions, it becomes possible to quantitatively investigate the 
similarities and differences between quantal and classical distribution functions in a broad class of smooth, conservative 
and strongly chaotic Hamiltonian systems. In particular, this paper displays (1) the rich transient behavior of the 
dynamics of distribution functions shared by quantum and classical dynamics before a QCC break time tf,, a quantity 
which scales with the logarithm of h, (2) a simple but enlightening description of the break regime of QCC in the time 
development of phase space structure, and (3) interesting QCC in low order statistical moments during a complex 
relaxation process with a time scale much larger than %. 

We note in passing that the distribution function strategy is also of considerable interest to the fundamental 
understanding of decoherence in quantum vs. classical mechanics. For example, the structure of the classical or 
quantum distribution functions determines properties of the system when coupled to an environment [T3.[l5| . Hence, 
QCC from a distribution function viewpoint is useful to the understanding of short-time quantum decoherence rates 
versus rates of classical entropy production Furthermore, QCC between quantal and classical distribution 

functions is closely related to the issue of decoherence in the presence of an environment that has a chaotic classical 
limit [HI- 

The model calculations in this paper are based on a strongly chaotic system: the quartic oscillator model . The 
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Hamiltonian is, in dimensionless scaled variables |19j], given by 

H( qi ,q 2 , PuP2 ) = f + f + ^hl + f (rf + rf). (1) 

When a = /3 or 3/3, this system is completely integrable. For very large values of a//3 the system is strongly chaotic, 
with the characteristic Lyapunov exponent being an order of magnitude larger than that of other conservative chaotic 
systems, e.g., the Hcnon-Heiles system. 

This paper is organized as follows. Section ITT1 briefly reviews the distribution function approach to QCC in chaotic 
systems. Special emphasis is put on a structure measure of classical and quantal distribution functions, and on our 
definitions of classical and quantal finite-time Lyapunov exponents. The paper then proceeds to present results in 
order of increasing time scales. In Sec. IIIII QCC is studied for initially positive-definite Wigner functions at early 
times. In Sec. IIVI a simple analytical understanding of the break time regime of QCC is provided, followed by 
supporting numerical results. We then consider, in Sec. QCC in some low-order statistical moments for much 
larger time scales. A brief summary is given in Sec. IVII 



II. LYAPUNOV EXPONENTS IN DISTRIBUTION FUNCTION DYNAMICS 



A. Classical distribution function dynamics 

Classical chaos is usually defined as the extreme sensitivity of trajectories to slight changes in initial conditions. 
Quantitatively, it is described by a positive maximum Lyapunov exponent or by the Kolomogorov entropy. Consider 
a conservative Hamiltonian system with two degrees of freedom with dimensionless canonical variables <?i, (?2,Pi,P2- 
A phase space point is characterized by a 4-dimensional column vector 7 = (91, 32>Pi>P2) • F° r brevity we introduce 

the antisymmetric matrix J = j J where and 1 are 2x2 zero and unit matrices, respectively. The canonical 
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equations of motion (i.e., Hamilton's equations) then read as 7 = JdH/dj. The sensitivity of classical trajectories 
to initial conditions is described by the stability matrix M,j [7(0), t] = dji(t)/djj(0), and the maximal Lyapunov 
exponent A is given by 

lim £ln(|M[7(0),tto|) = A[ 7 (0)], (2) 

where 77 is a vector in the tangent space. The stability matrix M is symplectic and its time evolution is governed by 
the differential equation 

M = J0M, (3) 

with 

M(0) = 1. (4) 

A simple application of Liouville's theorem leads to an alternative expression for the maximal Lyapunov exponent 
in terms of classical distribution function dynamics |lll Il2l . Let pt denote a well behaved classical probability 
distribution and £ denote an arbitrary infinitesimal vector in phase space. Then, from the incompressibility of the 
Liouvillc density function, one has 

pt[ 7 (i)] + Vpt[7(«)]M£ = po[j(0)} + Vp„[7(0)]£. (5) 

Further, using J 2 = —I and JM^J = — M _1 , one obtains 

MJVp [ 7 (0)]=JVp t [ 7 (i)]. (6) 

Substituting Eq. © into Eq. © yields 

A[ 7 (0),t] = lim iln|JVM7(*)]| 

t^oo t 

= hm iln|VM 7 (t)]|. (7) 
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Note that Eq. Q defines the Lyapunov exponent in terms of phase space distribution properties. Specifically, the 
faster the structure of p t increases, the larger the A. 

In accord with Ref. we define a particular ensemble average of the gradients of distribution functions as X2c 
i.e., 



X2c(t) 

Using Eqs. @) and JSJ, we have 



/|Vpt(7)| 2 rf7 
//»t(7)d7 . 



1/2 



J> t (7)V 2 p t ( 7 )d 7 



-,1/2 



I P 2 t{l)dl 



X2c(t) = 



J>o(7)d7 



|JM( 7s t)JVpo(7)l ^7 



1/2 



(8) 



(9) 



The physical significance of X2c(t) becomes apparent when one considers the Fourier transform of the distribution 
function. That is, suppose Pti'j) = [l/(27r) 4 ] / dkexp(27rik • j)p t (k), where pt(k) denotes the Fourier component 
evaluated at the 4-dimensional wave vector k. Then 



;^k 2 | Pt (k)i^ 
/dk| Pt (k)p 



(10) 



showing that \2c is the root-mean-square radius of the Fourier transform of the distribution function, and thus serves 
as a measure of classical phase space structure. That is, the larger the \2c, the more structured is pt( 7 ). 

For completely integrable dynamics there exists a special set of generalized coordinates: action variables Ix,Ii and 
angle variables #i, #2- In this representation, the Hamiltonian depends only on the action variables that are constants 
of motion. For such cases, Eq. (|3J) has the simple solution 



M=luS, 



(11) 



Substituting this explicit time dependence of the stability matrix into Eq. gives the following long-time behavior 

of X2c- 



lim 



X2c(t) 



t— >oo t 



I Plil)dl 



d 2 H 



dj 2 



JVp (7) 



1/2 



(12) 



Hence, for regular dynamics the structure of distribution functions, as measured by X2a asymptotically shows a linear 
time dependence in the action-angle representation. However, Xicif) may show a polynomial time dependence in 
other canonical representations. 

By contrast, for chaotic dynamics it was shown that |ll| 



lim - lnx2c(*) = A 2 , 

t— >oo t 



(13) 



where A2 is the so-called second order generalized maximal Lyapunov exponent. That is, in the chaotic case the root- 
mean-square Fourier radius of distribution functions increases, asymptotically, at an exponential rate of A2. Since 
a given resolution limit <5 corresponds to the inability to account for Fourier modes larger than 1/(5, chaos can be 
understood as a kind of exponential loss of accuracy, or of information, encoded in the Fourier basis expansion of the 
initial distribution function. 

By definition, the (generalized) Lyapunov exponent A2 is an asymptotic property, relevant as time goes to infinity. 
Realistically, however, it is finite-time properties of classical chaotic dynamics that are of real interest to the study 
of QCC. To this end it is useful to introduce a finite-time Lyapunov exponent. Based upon Eq. I|13|) . we define the 
finite-time Lyapunov exponent in terms of the average exponential increase rate of X2c(t) over time t: 



A 2c (t) 



1 



In 



X2c{t) 
X2c(0) 



(14) 



with 



lim X 2c (t) = A 2 . 

t — > + oo 



(15) 



We examine this quantity over various time scales later below. 
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B. Quantal analog of classical finite-time Lyapunov exponents 

From the ensemble point of view, QCC is best understood by comparing the classical Liouville equation with the 
quantum von Neumann equation in a phase space representation, e.g., the Wigner-Weyl representation. Specifically, 
given Eq. 10, which provides a quantitative diagnostic for characterizing classical chaos using classical distribution 
functions, it becomes straightforward to define the quantum analog of classical Lyapunov exponents using quantal 
distribution functions. By analogy with Eq. JHJ) we define the measure \2q f° r quantal phase space structure as 



X2q 



J\V P w h)\ 2 d 1 

/( P ^)2( 7 )d 7 



1/2 



(16) 



where p (7) is the Wigner function of a quantum state. Accordingly, in the Fourier space of quantal distribution 
function, 



a _ jdkk>r(k)p 

X2 « Jdk\ P f(kW 



(17) 



where pY {l) = [1/(2tt) 4 ] / dkexp(27rik • j)p^ (k). Hence, \2q is the root mean-square Fourier radius of the Wigner 
function. Further, finite-time Lyapunov exponents A2 9 (t) for quantum distribution function dynamics can be defined 
by direct analogy to \2c(t), i.e., 



A 29 (t) = -ln 



X2q(t) 
X2q(0) 



(18) 



Interestingly, \2q has an equivalent expression that is easier to handle. Suppose p is the density-matrix operator 
associated with the Wigner function p , and % is the operator associated with the classical canonical variable 7$ 
(e.g. qi,Pi). A simple calculation 0,0] then shows that 



xl 



, = 2 E 



Tr(p 2 jf - pjipji) 



h 2 Tr(p 2 ) 

Of particular interest is the pure state case, in which p 2 = p and as a result, 



x\ q 



■£«7?>-<7< 



(19) 



(20) 



where (•) represents ensemble expectation values. As shown later below, the analogous classical expression is far more 
complex. 

Consider then the quantum counterpart of Eq. (|15fl . Since Eq. <|19l) indicates that x\ q < 2 J2i Tr(p 2 'yi 2 )/?i 2 Trp 2 , 
we have that \2q has an upper bound for any bounded Hamiltonian system. As such, for fixed Ti and bounded systems, 



lim A 2 „(t) = 0. 

i— ^+00 



(21) 



This reproduces the widely accepted result that bounded quantum systems cannot exhibit chaos in the strict sense. 

However, as noted above, what is of interest to the study of QCC are the transient properties of \2q{t) versus \2 C {t)- 
Consider, for example, a two-degree-of-freedom system H = H{q\ 1 q2,pi,p2)- The quantum von Neumann equation 
in terms of the Wigner function p w is given by 



dp 



w 



Of 



{H,p w } 



E 

{h+h)>~L, odd 



(V2i) (/l+i2 ~ 1} d^+ l ^V( qi ,q 2 ) 8^ +l ^p W 



hiy. 



d<£d<£ 



(22) 



where the first term on the right hand side is the classical Poisson bracket, and the second term represents the sum 
over an infinite series of "quantum corrections". Consider the short time limit of X2 q (t). Using Eqs. (fTTHl and we 
have 



A 2g (0) 



1 



x E 

(h+l 2 )>l, odd 



d 2 H 



(^) J W>7 + J {Vp%Y d 7 
(h/2i)^+ 1 ^ d^+ l ^V{qu ffl) d^+^pg 



(23) 



where, obviously, the first term on the right hand side of Eq. (|23|) corresponds to the contribution from the classical 
Poisson bracket, and all other terms represent quantum "corrections". 

Some aspects of the distribution function strategy outlined above have been applied to the Arnold-cat map model 
|llL IT^ . where the stretching and contraction mechanism is uniform over the entire phase space. However, rich 
transient behavior of finite-time Lyapunov exponents, as implied in our derivation of the explicit state dependence of 
A2 9 (0), has not been explored thus far. The next section is devoted to both analytical and numerical studies on this 
subject. 



III. SHORT-TIME CORRESPONDENCE 



To examine the QCC in short-time dynamics it is useful to consider the classical analog of A2 g (0), i.e., A2 C (0) as 
the extreme short time limit of classical finite Lyapunov exponents. Using Eqs. ©, and JSJ, we have 



lim^ 
t^o dt 



(Vpo) T (S) J (Vpo)d7- 



Substituting Eq. $Q into Eq. (0 gives the zero time limit of \2 C (t) 



A 2c (0) = 



/(Vpo) T (0)J(Vpo)rf7 
J|Vp | 2 d7 



(24) 



(25) 



Here, A2 C (0) is seen to be the average of d 2 H/d^f 2 weighted by gradients of the initial distribution function. Not 
surprisingly, Eq. (|25|l resembles the first term on the right hand side of Eq. I|23|) . As A2 C (0) reflects an ensemble 
average of instantaneous density fluctuations, it depends strongly on the shape and location of the initial classical 
distribution function. Hence, as in the quantum case, there exists very rich transient behavior in the time development 
of phase space structure, an interesting feature that has often been ignored in previous QCC studies. 

To consider quantum effects induced solely by the dynamics suggests that we choose an initial quantal distribution 
function that is as classical as possible so that differences between classical and quantum dynamics evidently arise 
from the dynamics. This suggests that the initial Wigner function should be chosen as positive definite, so that it can 
be interpreted as a classical probability distribution. It is well known that for one-dimensional pure state dynamics 
the only positive-definite Wigner function is the Gaussian distribution function j^jj, El , which takes the following 
general form, 



o w — 



1 

— exp 
7m 



2^7% _, 2 



2r 



2 V 2 (l-r 2 ) + H(l-r 2 y/ 2iq " P) 



(26) 



where r, 77, q, p are parameters and where 

(q) = q, (p) = p, (q 2 - q 2 ) = »? 2 , 



(p 2 -p 2 ) 



4?7 2 (1 - r 2 ) 



(pq) - (p) (q) = 



hr 



2 (l- r 2)l/2' 



(27) 



This Gaussian form corresponds to the so called correlated coherent states \22\. whose coordinate representation is 
given by 



1 



(27T?7 2 ) 1 /4 



exp 



(l-r 2 )V2' 



-(a + a ) 

rj I 



where a is complex constant given by q/(2rf) + i [prj/h — qr/(2rjy/l — r 2 )] . In particular, for the case of 



(28) 



0, 



(q 2 -q 2 )(p 2 



V 



h /4, corresponding to the minimum- uncertainty-product state, i.e., the coherent state. 



For the two-degree-of-freedom system examined below one would choose two-dimensional Gaussian states. The 
initial quantal distribution function and classical distribution function po are thus chosen as the following, 

1 Pr 2 ,m,q->,P,- (29) 



Po 



Po 



Substituting this initial state into Eq. Q25J) an d approximating the average of the derivatives of ^(^1,52) as the 
derivative evaluated at the centroid of the Gaussian distribution, denoted d 2 V(q 1 ,q 2 )/^ f qh one S ets 
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It is seen that A2 C (0) depends strongly on both the shape parameters n, r 2 and the phase space locations of the initial 
distribution function. Two particular situations are worthy of note. First, for the minimum-uncertainty-product state 
(fi = r2 =0), i.e., the two-dimensional coherent state that is commonly used as initial states in QCC studies, Eq. 
(|3*0l gives A 2c (0) = 0. Second, when r, ; • (1 — d 2 V(q 1 ,q 2 )/dq 2 ) < for either i = 1 or i = 2, A2 C (0) can be negative. 
That is, for appropriate shape parameters r±, r-i and central coordinates qi,q 2 , the contraction mechanism associated 
with chaotic dynamics may initially dominate over the stretching mechanism, giving rise to a reduction in phase space 
structure. 

The quantum analog of this transient behavior can be examined by considering \2 q (t) in a similar fashion. Specifi- 
cally, substituting the initial positive-definte Wigner function l|29(l into Eq. (|23|l . keeping the leading order quantum 
correction term in the quantum Liouville equation, one obtains 



1 d 4 V 

A 2? (0) = A 2c (0)- 



167T f\VpW\2d 7 dqldql 



I dqid Pl \Vp™ uri ' 2 



' dq 2 dp 2 \Vp w 



(1— r 2 ) 1 / 2 J " rr 1 ,r 11 ,q 1 ,p 1 \ 

(31) 



(l — r 2)l/2 J yz L rr2,7i2,q 2 ,p 2 

Evidently, initial states with r% = r 2 = give A 2c (0) = A 2(3 (0) = 0. Thus, in this sense, the coherent state is the 
most classical state of the correlated coherent states; other types of initial states have a leading order quantum effect 
proportional to d 4 V '/ 'dq 2 dq 2 . After carrying out the integrals in Eq. i|31|) for given by Eq. 1)26(1. one sees that 

[A 2c (0) — A 2(? (0)] is proportional to the first power of h. Therefore, for relatively large h, depending upon the sign of 
r*i and r 2 , A 2(J (t) can be significantly larger or smaller than X 2c {t) at early times. In addition, if there is no quartic 
term in the potential, i.e., d 4 V/dq\dq2 = ( e -g-> m th° Henon-Heiles model), Eq. (|3~TJl shows that [A 2c (0) — A 2g (0)] 
is given by smaller terms that are proportional to higher powers of fi. 

Thus far we have only examined correspondence between A 2c (i) and A 2(7 (i) at t = 0. For nonzero times, one can 
utilize numerical methods to compare these quantities. For example, consider a coupled quartic oscillator system 
given by Eq. with a = 1.0 and (3 — 0.01, and for Ti = 0.5, 0.05, and 0.005. Classical calculations are done by 
Monte-Carlo methods based on Eq. 10, and quantum calculations use the FFT split operator technique 23]. For 
each case we examine three sets of initial distributions, i.e., (a) r\ = r 2 = 0, (b) n = r 2 = 0.6, and (c) r\ = r 2 = —0.6. 
For all three cases rji = r)2 = ^/h/2, and the centroid of the initial state is fixed at q 1 = 0.40, q 2 — 0.60, p 1 = 0.50, 
p 2 = 0.414. Note that this initial location of the Gaussian distribution gives (1 — d 2 V/dq\) > and (1 — d 2 V/dq 2 ) > 0, 
a fact that is relevant to the discussions below. 

Figures Hl3l display comparisons between X2c{t) and X2q(t)- Consider Fig. for the case of Ti = 0.5. Here 
d 4 Vjdq\&q\ — 2.0, and one expects that A 2g (0) can deviate considerably from A 2c (0) for nonzero r\ and r 2 . Further, 
the expectation [from Eq. I|31|) ] is that the initial quantum correction should be positive when T\, r 2 are negative, and 
negative when ri,r 2 are positive. All these predictions are confirmed nicely in our computations. In particular: (i) in 
Fig. ^ ( r i — r 2 = 0), both X2c(t) and X2q(t) assume an initial zero slope; (ii) in Fig. ^ ( r i = r 2 = 0.6), the positive 
slope of X2q(t) at t = is seen to be smaller than that of X2c(t); and (iii) in Fig. ^ ( r i = r 2 = —0.6), it is indeed 
seen that at very short times X2q(t) > X2c(t), both of which have negative initial slopes. For all three situations, the 
quantum effects are so large that QCC is essentially lost at very short times. 

Figure shows the corresponding results (with different abscissa scale) with the value of h decreased by a factor 
of 10. This case is different from that in Fig. ^ insofar as excellent QCC is seen for short times for all three initial 
conditions. In particular, in Fig. [2b the time dependence of both X2c(t) and X2q(t) is seen to be very close to 
exponential for t < 1.5. However, this is partially because the initial slope of the curves [i.e., A 2g (0) and A 2c (0)] in Fig. 
|2t> both happen to be close to the average slope (i.e., finite-time Lyapunov exponents) for longer times (e.g., t ~ 4), 
In Fig. [2j the transient behavior of the finite-time Lyapunov exponents assumes a completely different nature, i.e., 
the increase of both X2c(t) and X2q{t) is significantly suppressed for times up to t = 0.5. Variations in the finite-time 
Lyapunov exponents can be further seen by comparing X2c(t = 4.0) and X2q{t = 4.0) in Fig. 03 to those in Fig. [5^, 
and Fig. they can differ by a factor as large as 1.5. Note also that Fig. [5] shows that the break regime of QCC 
(i.e., the time when classical and quantum dynamics no longer agree) is between t — 1.0 and t = 1.5, which is on 
the order of one average period of motion. With increasing time, the agreement between X2q(t) and X2c(t) in Fig. [21 
worsens: the classical phase space structure is seen to increase exponentially on the average, whereas there is no clear 
sign of a similar exponential increase in the quantum distribution dynamics. 

As an aside we note that since the structure of quantal distribution functions may determine the early-time deco- 
herence rate if the quantum system is open |14| . these variations in finite-time Lyapunov exponents imply that even 
in classically strongly chaotic systems it is still possible to pick out some initial states which are relatively robust to 
decoherence effects. 

As shown in Fig. |3 with a further large decrease of h, the break regime of QCC is considerably shifted, to 
t = 2.0 — 3.0. The transient behavior caused by different initial states still persists for t < 0.5. After the transient 
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'0 0.25 0.5 0.75 1 

time 




3.1 - 



3 0.25 0.5 0.75 1 

time 

FIG. 1: Time dependence of X2c(i) (dashed line) and X2q{t) (solid line) for three initial Gaussian distributions (h — 0.5). The 
three panels correspond to (a) ri = ri — 0, (b) ri = ri = 0.6, (c) ri = V2 = —0.6. Note that the initial slope of these solid and 
dashed lines are zero, positive, and negative in (a), (b) and (c), respectively. All variables are in dimensionless units. 



period and before the break time, all quantal distribution functions in Figs. [2K _C emanating from the three different 
initial conditions are seen to undergo an exponential increase of structure on the average, in excellent agreement with 
the behavior of classical distribution functions. 
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FIG. 2: Same as Fig. ^except fi = 0.05 and that both X2c(t) and Xtq(t) are plotted in the logarithmic scale. The average slope 
of the curves (up to a certain time) is indicative of the finite-time Lyapunov exponents \2c{t) and \2 q {t). The break regime of 
QCC is clearly between t — 1.0 and t = 1.5. All variables are in dimensionless units. 



IV. THE BREAK REGIME FOR CORRESPONDENCE 



Understanding the QCC break regime is a central problem in the study of correspondence. From the viewpoint 
adopted here, i.e., via phase space distribution functions, the break regime is when the classical distribution function 
begins to develop significantly different structure from the quantal distribution function, i.e., X2c{t) begins to deviate 
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FIG. 3: Same as Fig. |U except h = 0.005. The break regime of QCC is clearly between t = 2.0 and t = 2.5. For 0.5 < t < 2.5 
and for all three initial conditions in (a), (b), and (c), there is an exponential increase of phase space structure for both classical 
and quantum dynamics. All variables are in dimensionless units. 



appreciably from X2q(t)- 

It is tempting to try to use the quantum Liouville equation [Eq. I)22[l] to study the breakdown of QCC with time. 
However, this is not an easy task, since, for nonzero times, one cannot analytically deal with distribution function 
dynamics. Rather, we present a simple description of the QCC break regime by first re-interpreting xicW an d then 
comparing it with xi q (t). As will be quite clear, our description also relates the distribution function strategy to the 
trajectory viewpoint of chaos. 
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It is convenient to restrict the discussion to the case of the minimum-uncertainty-product state, although the 
following considerations can be readily extended to the most general case. Further, we focus on pure state dynamics 
because mixed state dynamics simply makes the study of QCC more complicated, without adding new physical insight. 

We begin with the Taylor series expansion of the stability matrix M[7(0),i] around 7(0); i.e., 

M kim ,t] = M tt [7(0),*]+E^(0) y^'* ] 



m mn 

+O((M0)) 3 ), (32) 

where we have defined A k i, Bkim, and Gkimn, as the zeroth, first, and second order derivatives of the stability matrix 
at 7(0), respectively. The summation indices m,n, etc. run from 1 to 4 for a two-degree-of- freedom system, and 
£7(0) = [7(0) — 7(0)]. The time dependence of Aki, B k i m , and Ckimn is governed by the following set of first order 
differential equations together with the canonical equations for classical trajectories: 

dA kl ^ T d 2 H 



dt jjj dji'd^k' 
dB klm T /v- d 3 H A A , d 2 H 



j. / " HI \ / . f~\ r\ i") 71 ill K I ' r\ r\ 




B k 



k'lri 



d Ckimn _ \ - T d_ H d H 

,. / 'Jkl' I / , o r\ <~x ^m'n&k'lm "T ^ ^k'lnin 

dt frf dji>dj k iOJ m > oil' oik' 

d 3 H 



k'l' n'rn' 



dji'djk'di 
& 

djl'djk'dlrn'djn 



(B 



] JwT, o o o Ak'iA n ' m A m i n . (33) 



Consider two initial points in the phase space, 7(0) and 7(0). They generate two classical trajectories, denoted by 
7(t) and j(t), respectively. From the trajectory viewpoint, of most interest is the difference Sjj(t) = [jj(t) — Jj(t)], 
which is a differentiable function of 7(0), Sj(0), and t. Fixing 7(0), and thus j(t), gives a reference trajec- 
tory. Consider now the Taylor series expansion of Sjj(t) in terms of 57(0). Obviously Sjj(t) — if S-f(0) = 0, 
d5jj(t)/d6'yi(0)\s~/(o)=o = Aji, d 2 5jj(t)/d5jk(Q)d5ji(0)\s'y(o)=o = Bjki , and so on. In fact, the nth order derivative 
of S~fj(t) with respect to S-f(0) is simply given by the (n — l)th order derivative of the stability matrix M[t"(0),£] with 
respect to 7(0). Based on this observation, one obtains the following relation: 



5^*71(0) = S 7j (t)-^S lk (0)6 7l (0)B jkl 

Ik 



~\ E ^k(0)S 7l (0)S 7m (0)C jkln - O ((5 7 (0)) 4 ) • (34) 

jklm 

Consider now the initial classical distribution function corresponding to the coherent state 



: ' ,0, - T,l '' ,] = U.V la a ' ' M> 



(gi ~ 9i) 2 _ (92 ~ q\) 
2 < 2< 



{pi-pi) 2 (P2-P2) 2 ' ^ 



2 < 2 < 



After lengthy calculations (see Appendix) using Eqs. (JSJ, (|3~4"|) and (J2U, we obtain an enlightening expression 
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for xlc(«): 



xlc(*) = |rE[«^(*)) 2 >°- (7i(*)>§] 

E Qs^B^ + \A jh C jmi ) («7k(0)*Tv(0)*7l(0)*7J'(0)>o 
ffefc'W ^ 'J 



2 



jkk'W v 

+ io((5 7 (0)) e 
7i 



(36) 



Here (-)o represents the average over the initial classical ensemble. Note that the factor h in the above equation 
is due to the fact that a qi a pi — cr q2 a P2 — Ti/2 for any classical distribution function corresponding to a minimum- 
uncertainty-product state. 

Equation (|36l) provides a quantitative connection between the development of phase space structure and the insta- 
bility characteristics of classical trajectories. Further, it allows for a closer examination of when and how classical 
distribution functions begin to develop structure that differs from quantal distribution functions. Consider then, for 
simplicity, a special case in which the initial coordinate variances are identical to the initial momentum variances, i.e., 

(Tq, — (J 0l — <j q 



<J P2 . One obtains 



xUt) = T2 E [<(7,W) 2 )o- (7i(*))g] + f(t)+0(h), 



(37) 



where f(t) is given by 
/(*) 



/l 2 
E ( ^BjkiBjk'V + -^AjkCjyw 



jkk' IV 

3 1 

-^Skk'w + j [(Skk'dw + Ski'Sik')(l - $ki) + Skifa>i'(l - fak-)] 

E JkmJk'm' [BjklBjk'V + Aj k Cjk>U>) 
jkk' 11' mm' 

3 1 

-^mm'W + ^ [(fimm'Sll' + 5 m i'Si m ')(l — S m i) + S m iS m 'i'(l — S mm ')] 



(38) 



This expression affords new perspectives into QCC. Recall first, as shown in the previous section [see Eq. 1)20(1 ]. that 
for pure state dynamics the quantal measure \2q (0 °f phase space structure can be expressed in terms of the sum of 
second order moments. This is more or less a consequence of quantization: the smallest scale in coordinates is related 
to the largest momentum component of the wavefunction, and vice versa. Clearly, this is in general not the case for 
classical distribution functions. Nevertheless, Eqs. I|36|) an d l|37|l indicate that, to lowest order, X^ c {t) can still be 
related to the sum of second order statistical moments. Specifically, as shown in Eq. I|37(l . for an initially symmetric 
coherent state the first term contributing to xLW ^ s gi ven by the sum of second order moments divided by ?i 2 /2, 
which is exactly the same as the result in Eq. 12U|) for x| g (t). In addition, xicW contains additional contributions 
absent in the quantum dynamics. For example, the leading order correction term f(t) is independent ofh. As shown 
by Eq. (|38|l . this term is determined by the stability characteristics A/fc, Bjk m , and Cjkmn associated with the 
trajectory starting from the centroid of the initial Gaussian distribution. These results show intriguing similarities 
and differences between X2c(t) and X2q(t). 

As a simple example of Eq. i|37|) i consider a quadratic Hamiltonian systems such as the harmonic oscillator or the 
inverted harmonic oscillator system, where classical and quantum mechanics are expected to agree. Here the time 
evolution is a linear canonical transformation in phase space; hence the stability matrix elements M.jk[j(0), t] do not 
depend on 7(0) and Bjki, Cjkim and all other higher order derivatives of the stability matrix with respect to 7(0) are 
zero. Thus, in such linear systems, xicW is precisely given by the sum of some second order moments divided by 
h 2 /2, in perfect correspondence with xlqW [ see Eq. lj2"U|) ]. 

For the case of chaotic systems, one can estimate that the stability matrix increases exponentially, i.e., A^i ~ 
exp(At), with A being the average exponential increase rate up to time t. Likewise, one expects B^im ~ exp(At), 
Ckimn ~ exp(Ai), etc. Equation then suggests that f it) ~ /oexp(2A£). By contrast, the first term on the right 
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hand side of Eq. (|37|l is determined by second-order moments and cannot increase for all time for bounded systems. 
Indeed, assuming that the characteristic magnitude of the second order moments [((7j(i)) 2 )o ~ (7i(*))o] U = 1> 2, 3, 4) 
is given by O 2 , then the first term on the right hand side of Eq. 1)3 7|) would be bounded by 80 2 /?i 2 . Thus, in Eq. 
(|37|l the f{t) term will be comparable to the preceding term after a time tb, approximately given by 



t b = -In 

A comparison between Eqs. I|2(J|) and l|37|l suggests that tb corresponds to the time scale after which classical descrip- 
tions of phase space structure no longer agree with quantum results. Thus, tb can be identified as a logarithmic break 
time of QCC. This result is consistent with previous studies on the QCC break time using different approaches [24| . 
Note that, since Eq. I|39|) involves classical variables only, one can calculate tb without the need for any quantum 
calculations. 

The origin of the classical-quantum difference lies in the f(t) term in Eq. 1)37(1. This term does not have a quantal 
analog since it reflects classical phase space structure that is beyond the resolution limit of quantal distribution 
functions [2^|. To demonstrate the role of f(t) we compare [xlcW — xlgW] from a direct calculation to f(t) given by 
Eq. 1|38[) . Specifically, f(t) is obtained from Eq. (|38|l by numerically computing the instability characteristics Ajk, 
Bjki and Cjkim via Eq. (|33|l . Two cases with differing values of U have been examined. Results are shown in Figs. 
0Ji and^J) that display the comparison between [xf c (t) — xi q {t)} and f(t) for both cases. The agreement is excellent, 
confirming the role of f(t) in determining the break regime of QCC. 



(39) 




time time 



FIG. 4: Time dependence of (xlc — xi 9 ) compared with f(t) in our analytical considerations [see Eqs. 13 7H and 13811 in the 
text]. The initial condition is the minimum-uncertainty-product state considered in Sec. Dashed lines denote (xic — xi 9 ) 
based on direct calculations in Sec. IIIII and solid lines denote f(t) obtained by directly integrating Eq. I|38^ . Ti equals 0.05 in 
(a) and equals 0.005 in (b). The agreement between (xic — X2q) and f(t) for the QCC break regime is excellent. All variables 
are in dimensionless units. 

Figure 0] also shows some tiny discrepancies between [xi c (*) — X2 g M] an d /(*) (especially at later times), corre- 
sponding to higher order contributions in Eq. I|37|) . Nevertheless, the physics is still the same; i.e., the sum of f(t) and 
these higher order terms measures very fine phase space structure that is unresolvable by the quantum distribution 
function. 



V. LONG AFTER THE BREAK REGIME 



The logarithmic break time tb of QCC introduces many issues regarding the relationship between quantum dynamics 
and classical mechanics. For instance, Zurek et al. |26| showed that a logarithmic break time (somewhat different 
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from that obtained above) can be counter-intuitively short even for macroscopic objects. This being the case, our 
everyday experience in a macroscopic classical world full of chaotic events seems incompatible with the notion that 
classical physics is a large-quantum-number limit of quantum mechanics. Likewise, the smallness of the QCC break 
time seems to imply that classical physics cannot play a role in nonlinear molecular dynamics, contradicting the fact 
that classical physics often works well in many dynamics simulations. To at least partially resolve these puzzles, this 
section attemps to explore the implications of the quantum-classical discrepancy in phase space structure for ensemble 
statistics. 

As already implied by a comparison between Eqs. I|20[l and H37|). what is directly responsible for the logarithmic 
break time is not quantum classical difference in expectation values, but simply the exponential increase of f(t) which 
reflects the richness of fragmentation of classical distribution functions. As such, it is interesting to examine QCC in 
terms of some observables. Figure |S] displays the time dependence of four variances, i.e., (qf) — (q\) 2 , (q%) — ((fe) 2 , 
(Pi) ~ (Pi) 2 > an d {p 2 ) — (P2) 2 , for both classical and quantum dynamics. The initial state corresponds to that used in 
Fig. and Fig. ^J) and h still equals 0.005, allowing direct comparison with the previous results. 



1.8 1 1 1 1 1 1 1.2 




time time 



0.2 




time 



0.2 




time 



FIG. 5: Quantum classical comparison of the time dependence of four second-order statistical moments, i.e., {q 2 ) — (qi) 2 in 
( a )i (<&) ~ (Q2} 2 in (b), (pi} — (pi) 2 in (c), and {p 2 ) — (P2} 2 in (d), for times much larger than the logarithmic break time tb- 
The initial distribution function corresponds to that in Figs. and^Ja, and Ti = 0.005. Dashed lines denote classical results, 
solid lines are quantum results. The QCC shown here during the complex relaxation process (0 < t < 12.0) is remarkable. All 
variables are in dimensionless units. 

Interestingly, results in Fig. [S]show that QCC in these second order moments is excellent for times up to t « 12.0, 
which is much longer than the break time tf, ~ 3.0 identified in Fig. |3| In particular, for times less than t = 2.5, both 
quantal and classical variances are seen to increase very rapidly, in exactly the same manner. This rapid increase is 
exponential in nature because, in this regime, the sum of these variances (divided by 7i 2 /2) give x 2 c ( or x!g)> which 
indeed increases exponentially on the average (see Fig. [3^). This confirms a published suggestion that the initial rapid 
increase of both quantal and classical variances could be used to identify classical chaos from quantum dynamics Q • 
For times 2.5 < t < 12.0 during which excellent QCC persists, the second order moments display a complex evolution. 
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For example, there is clearly a diffusive regime between t = 5.0 and t = 7.0 with very large diffusion constants. 
After that, several very strong oscillations (particularly in Figs. 0>d) can be observed, suggesting that both the 
quantal and classical distribution functions alternate between a high degree of derealization and a certain degree of 
localization over the entire accessible phase space (this oscillatory behavior also suggests that the QCC time scale 
here should scale as a power of h, rather than scale as InTi). The impressive QCC in this complex relaxation pattern 
is in sharp contrast to the results in Figs. and^J), where quantum classical differences, as quantitatively described 
by [x!e(£) — x!g(^)] an d predicted by f(t), are already very large for 3.5 < t < 4.0. For even later times (t > 12.0), Fig. 
[3] shows that quantum classical discrepancies begin to show up, but still with very similar trends in the oscillations. 

These results indicate that the logarithmic break time may be quite irrelevant to QCC when viewed from the 
perspective of some low order statistical moments |2^. This is the case since considering low order statistical moments 
is equivalent to projecting the quantal and classical distribution functions onto a much smaller subspace. In this 
projection, or coarse-graining procedure, all the information encoded in an infinite number of higher order statistical 
moments is lost. An accurate QCC is thus restored due to the loss of detailed information. 

The results here also give further support to the idea of applying classical propagation methods to quantum 
distribution functions, as a means of approximating the true quantum mechanics |27j. That is, chaos may not cause 
a rapid failure of classical dynamics simulations if one is only interested in expectation values of some observables. 
More importantly, the results provide more insights into a fundamental issue regarding the role of decoherence in 
QCC [29l l30| . On the one hand, our results here show that during a relaxation process whose time scale is much 
larger than the logarithmic break time, decoherence may still be unnecessary in order to ensure QCC in low order 
statistical moments. Thus, at least in our model system, decoherence effects on correspondence in expectation values 
of low order statistical moments, if any, should be considered after the relaxation process is essentially completed, 
e.g., after t = 12.0 in Fig. [SJ This point supports the argument of Casati et al. |3C|. extending their considerations 
from one-dimensional kicked systems to conservative systems, and agrees with the recent work kl by Emerson and 
Ballentine. On the other hand, as shown in our previous work using the same model system |l9j |. QCC for much 
larger time scales can be much improved with the introduction of decoherence. 



VI. SUMMARY 



We have studied the issue of QCC in conservative chaotic systems in detail using a phase space distribution function 
approach. The nature of QCC in three different regimes is exposed. In particular: (i) the short time increase rate 
of phase space structure is studied in connection with finite-time Lyapunov exponents; (ii) a simple and enlightening 
description of the break regime of QCC is obtained, by successfully accounting for the classical phase space structure 
which is beyond the quantum description; and (iii) excellent QCC in low order statistical moments is shown to persist 
during a complex relaxation process, whose time scale is much larger than the logarithmic break time. 
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APPENDIX A: DERIVATION OF EQ. (36) 



In this appendix we give a rather detailed derivation of Eq. I)30H using Eqs. JjJJl, l|32|) . i|34|) and (|35|l . We begin 
with the derivation by substituting Eq. (|3"2|> and the initial distribution function p [y(Q), 7(0)] [see Eq. |3"5)l] into Eq. 
(jUJ, evaluating the first derivatives of po [7(0), 7(0)], and finally rescaling all integration variables by a factor of l/\/2. 
One then obtains 



xl(t) 



d 7 (0W7(0),7(0)] 



E 



E 



+ E J C jkmnSl m (0)S ln (0) + O ((<5 7 (0)) 3 



(Al) 
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where 



Oil 



«3 



Pi(0)-Pi(0) 



in 



a 2 



9i(0)-9i(0) 



P2 (0)-p 2 (0) 
«fe(0)-g 2 (0) 



(A2) 



Since A,fc, Bjkrm and Cjkmn do not depend upon 7(0), and only the even functions of (57(0) will contribute when 
integrated over po [7(0), 7(0)], Eq. (|A1|) can be further reduced to 

'mm' + -BjkmBjk'm* J {(^kO-k' Sy m (0)Sj m * (0)) 

j fcfc' mm' 



~ J d 7 po[7(0),7(0]) (EIE A ^| 2 ) +O(a 2 (M0)) 4 ), 



(A3) 



where (-)o means the ensemble average over the initial Gaussian distribution. 

To further simplify the expression for xlcM' we make a change of the integration variables in the second term on 
the right side of Eq. I|A3|) . i.e., 



(Pi -Pi) 



-(q-Qi), -(P2-P2) 



-(32 - 3 2 )i 



'</2 



(qi-Qi) -> — (P1-P1), (92 - g 2 ) — (Pa -p 2 )- 



(A4) 



Note that the Jacobi matrix of this transformation is a unity matrix and it has no effect on the form of po [7(0), 7(0)]. 
As a result of this coordinate transformation we have otf. — > 8^(0)/ (cr qi <T pi ) =5^k{Q) / (a q2 a P2 ). This makes it possible 
to reexpress the second term on the right side of Eq. (|A3|) by use of Eq. (|34J) . Specifically, 



jkOLk\ 



\ f d 7 (0)p o [ 7 (0), 7 (0)] £lE^ 

V 3 k 

d 7 (0)p [ 7 (0),7(0)] (ElE A ^( )| S 



2 K 1 CTp 1 ) 2 
1 



E «(7;(*)) 2 >o - (7i(*)>o) - £^W)*7*(0)*yi(0)> 



2 0gi°>i) 2 

+ 7 E ^fei^fe'/'(' 5 7fc(0)^ 7 fc'(0) ( 57 / (0)5 7 p(0)) 



jklk' V 



\ £ a j - Wro <« 7i (t)J 7Jk (Q)57 l (0)*7«.(0)>o + O ((MO)) 6 ) 



jklm 



(A5) 



The second and the fourth terms in the above expression arc a linear function of Sj(t) and they can be further 
transformed into some functions of <5 7 (0) by using Eq. ffi^i a second time. We then get 

\ J d 7 (ow7(o), 7 (o)] I EiE^^i 

\ 3 

■E[<^(*)) 8 )o-(7i(t)>?] 

2 (jB^S^ + ^j*C ifc /„/)(*7 fc (0)*7 fc i(0)*7i(0)*>(0))o 
rO((<5 7 (0)) 6 ). 



2(<j qi a Pl ) 



2(ct 0i ct„,) 2 ^ M~ 



(^l^pi) 2 



(A6) 
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Finally, inserting Eq. HA6f) into Eq. (|A3|) yields Eq. 
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